Conformational Changes and Unfolding of ß-Amyloid Substrates in the Active Site of γ-Secretase.

Alzheimer's disease (AD) is the leading cause of dementia and is characterized by a presence of amyloid plaques, composed mostly of the amyloid-ß (Aß) peptides, in the brains of AD patients. The peptides are generated from the amyloid precursor protein (APP), which undergoes a sequence of cleavages, referred as trimming, performed by γ-secretase. Here, we investigated conformational changes in a series of ß-amyloid substrates (from less and more amyloidogenic pathways) in the active site of presenilin-1, the catalytic subunit of γ-secretase. The substrates are trimmed every three residues, finally leading to Aß40 and Aß42, which are the major components of amyloid plaques. To study conformational changes, we employed all-atom molecular dynamics simulations, while for unfolding, we used steered molecular dynamics simulations in an implicit membrane-water environment to accelerate changes. We have found substantial differences in the flexibility of extended C-terminal parts between more and less amyloidogenic pathway substrates. We also propose that the positively charged residues of presenilin-1 may facilitate the stretching and unfolding of substrates. The calculated forces and work/energy of pulling were exceptionally high for Aß40, indicating why trimming of this substrate is so infrequent.

Interaction of ROMK2 channel with lipid kinases DGKE and AGK: Potential channel activation by localized anionic lipid synthesis.

In this study, we utilized enzyme-catalyzed proximity labeling with the engineered promiscuous biotin ligase Turbo-ID to identify the proxisome of the ROMK2 channel. This channel resides in various cellular membrane compartments of the cell including the plasma membrane, endoplasmic reticulum and mitochondria. Within mitochondria, ROMK2 has been suggested as a pore-forming subunit of mitochondrial ATP-regulated potassium channel (mitoKATP). We found that ROMK2 proxisome in addition to previously known protein partners included two lipid kinases: acylglycerol kinase (AGK) and diacylglycerol kinase ε (DGKE), which are localized in mitochondria and the endoplasmic reticulum, respectively. Through co-immunoprecipitation, we confirmed that these two kinases are present in complexes with ROMK2 channels. Additionally, we found that the products of AGK and DGKE, lysophosphatidic acid (LPA) and phosphatidic acid (PA), stimulated the activity of ROMK2 channels in artificial lipid bilayers. Our molecular docking studies revealed the presence of acidic lipid binding sites in the ROMK2 channel, similar to those previously identified in Kir2 channels. Based on these findings, we propose a model wherein localized lipid synthesis, mediated by channel-bound lipid kinases, contributes to the regulation of ROMK2 activity within distinct intracellular compartments, such as mitochondria and the endoplasmic reticulum.

Investigation of Peptides for Molecular Recognition of C-Reactive Protein-Theoretical and Experimental Studies.

We investigate the interactions between C-reactive protein (CRP) and new CRP-binding peptide materials using experimental (biological and physicochemical) methods with the support of theoretical simulations (computational modeling analysis). Three specific CRP-binding peptides (P2, P3, and P9) derived from an M13 bacteriophage have been identified using phage-display technology. The binding efficiency of the peptides exposed on phages toward the CRP protein was demonstrated via biological methods. Fibers of the selected phages/peptides interact differently due to different compositions of amino acid sequences on the exposed peptides, which was confirmed by transmission electron microscopy. Numerical and experimental studies consistently showed that the P3 peptide is the best CRP binder. A combination of theoretical and experimental methods demonstrates that identifying the best binder can be performed simply, cheaply, and fast. Such an approach has not been reported previously for peptide screening and demonstrates a new trend in science where calculations can replace or support laborious experimental techniques. Finally, the best CRP binder─the P3 peptide─was used for CRP recognition on silicate-modified indium tin oxide-coated glass electrodes. The obtained electrodes exhibit a wide range of operation (1.0-100 µg mL-1) with a detection limit (LOD = 3σ/S) of 0.34 µg mL-1. Moreover, the dissociation constant Kd of 4.2 ± 0.144 µg mL-1 (35 ± 1.2 nM) was evaluated from the change in the current. The selectivity of the obtained electrode was demonstrated in the presence of three interfering proteins. These results prove that the presented P3 peptide is a potential candidate as a receptor for CRP, which can replace specific antibodies.

GS-SMD server for steered molecular dynamics of peptide substrates in the active site of the γ-secretase complex.

Despite recent advances in research, the mechanism of Alzheimer's disease is not fully understood yet. Understanding the process of cleavage and then trimming of peptide substrates, can help selectively block γ-secretase (GS) to stop overproduction of the amyloidogenic products. Our GS-SMD server (https://gs-smd.biomodellab.eu/) allows cleaving and unfolding of all currently known GS substrates (more than 170 peptide substrates). The substrate structure is obtained by threading of the substrate sequence into the known structure of GS complex. The simulations are performed in an implicit water-membrane environment so they are performed rather quickly, 2-6 h per job, depending on the mode of calculations (part of GS complex or the whole structure). It is also possible to introduce mutations to the substrate and GS and pull any part of the substrate in any direction using the steered molecular dynamics (SMD) simulations with constant velocity. The obtained trajectories are visualized and analyzed in the interactive way. One can also compare multiple simulations using the interaction frequency analysis. GS-SMD server can be useful for revealing mechanisms of substrate unfolding and role of mutations in this process.

Protein Homology Modeling for Effective Drug Design.

The effective drug design, especially for combating the multi-drug-resistant bacterial pathogens, requires more and more sophisticated procedures to obtain novel lead-like compounds. New classes of enzymes should be explored, especially those that help bacteria overcome existing treatments. The homology modeling is useful in obtaining the models of new enzymes; however, the active sites of them are sometimes present in closed conformations in the crystal structures, not suitable for drug design purposes. In such difficult cases, the combination of homology modeling, molecular dynamics simulations, and fragment screening can give satisfactory results.

Specificities of Protein Homology Modeling for Allosteric Drug Design.

The allosteric binding sites are usually located in the flexible areas of proteins, which are hardly visible in the crystal structures. However, there are notable exceptions like allosteric sites in receptors in class B and C of GPCRs, which are located within a well-defined bundle of transmembrane helices. Class B and C evolved from class A and even after swapping of orthosteric and allosteric sites the central binding site persisted and it can be used for easy design of allosteric drugs. However, studying the ligand binding to the allosteric sites in the most populated class A of GPCRs is still a challenge, since they are located mostly in unresolved parts of the receptor's structure, and especially N-terminus. This chapter provides an example of cannabinoid CB1 receptor N-terminal homology modeling, ligand-guided modeling of the allosteric site in GABA receptor, as well as C-linker modeling in the potassium ion channels where the allosteric phospholipid ligand PIP2 is bound.

Antioxidant activity of novel diosgenin derivatives: Synthesis, biological evaluation, and in silico ADME prediction.

A series of novel diosgenin (DSG) derivatives has been synthesized and tested in vitro for their antioxidant activity. Initially, four analogues have been evaluated for their cytotoxicity using normal human skin fibroblast (NHDF) as model cells. As a result, 84% of NHDF cells were still alive at 5 µM, so these compounds can be considered as innoxious to fibroblasts at this concentration. Then, hemolytic activity against human erythrocytes was studied in order to evaluate the potential impact of tested compounds against normal host cells. The result < 5% of hemolysis rates suggest no lytic activity for most compounds. After that, the main test - evaluation the antioxidant effect of DSG and its new derivatives against lipid peroxidation in the o/w emulsion model - was performed. The most promising compound (8) exhibited the significant antioxidant activity and the biocompatibility towards normal human dermal fibroblasts and red bloods cells. This p-aminobenzoic derivative revealed 61.6% blocking of induced lipid oxidation. Furthermore, eleven predicted ADME properties were predicted for all tested compounds and revealed that they are in compliance with drug-likeness criteria.

COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments.

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.

Evaluation of immune evasion in SARS-CoV-2 Delta and Omicron variants.

Emerging SARS-CoV-2 variants with higher transmissibility and immune escape remain a persistent threat across the globe. This is evident from the recent outbreaks of the Delta (B.1.617.2) and Omicron variants. These variants have originated from different continents and spread across the globe. In this study, we explored the genomic and structural basis of these variants for their lineage defining mutations of the spike protein through computational analysis, protein modeling, and molecular dynamic (MD) simulations. We further experimentally validated the importance of these deletion mutants for their immune escape using a pseudovirus-based neutralization assay, and an antibody (4A8) that binds directly to the spike protein's NTD. Delta variant with the deletion and mutations in the NTD revealed a better rigidity and reduced flexibility as compared to the wild-type spike protein (Wuhan isolate). Furthermore, computational studies of 4A8 monoclonal antibody (mAb) revealed a reduced binding of Delta variant compared to the wild-type strain. Similarly, the MD simulation data and virus neutralization assays revealed that the Omicron also exhibits immune escape, as antigenic beta-sheets appear to be disrupted. The results of the present study demonstrate the higher possibility of immune escape and thereby achieved better fitness advantages by the Delta and Omicron variants, which warrants further demonstrations through experimental evidences. Our study, based on in-silico computational modelling, simulations, and pseudovirus-based neutralization assay, highlighted and identified the probable mechanism through which the Delta and Omicron variants are more pathogenically evolved with higher transmissibility as compared to the wild-type strain.

Flavonoid quercetin abolish paxilline inhibition of the mitochondrial BKCa channel.

Large-conductance calcium-regulated potassium channel (BKCa) is known to play an important role in physiological and pathological processes. Despite the BKCa channel being encoded by one gene, this channel has been found to be located not only in the cell membrane but also in the membranes of intracellular compartments, such as in the inner mitochondrial membrane. With some differences, the mitochondrial BKCa (mitoBKCa) channel has been shown to be activated or inhibited by both synthetic and natural compounds. One of them, paxilline, has been considered to be a canonical blocker of this channel. In the previous study, we showed that the natural origin substance quercetin activates the mitoBKCa channel at ten times lower the concentration compared to channel present in the plasma membrane. Here, using the patch-clamp technique, we report that after inhibition of mitoBKCa channels by paxilline, quercetin activates these channels, indicating a paxilline and quercetin binding competition in the regulation of the mitoBKCa channel. To support our hypothesis, we used an analog of quercetin - isorhamnetin, a substance with one substituent changed. Isorhamnetin has no effect on the mitoBKCa channel activity, and after its application, paxilline fully inhibits the channel. Additionally, the molecular modeling studies were used. The results of docking quercetin and paxilline to the BKCa channel suggest that paxilline cannot bind after activation of the channel with quercetin. It seems that the likely mechanism of this phenomenon is the formation of spatial hindrance by quercetin. The results obtained shed a completely new, groundbreaking in the paxilline context, light on the current knowledge about mitochondrial potassium channel regulation.

SOD1 mutations associated with amyotrophic lateral sclerosis analysis of variant severity.

Mutations in superoxide dismutase 1 gene (SOD1) are linked to amyotrophic lateral sclerosis (ALS), a neurodegenerative disorder predominantly affecting upper and lower motor neurons. The clinical phenotype of ALS shows inter- and intrafamilial heterogeneity. The aim of the study was to analyze the relations between individual SOD1 mutations and the clinical presentation using in silico methods to assess the SOD1 mutations severity. We identified SOD1 causative variants in a group of 915 prospectively tested consecutive Polish ALS patients from a neuromuscular clinical center, performed molecular modeling of mutated SOD1 proteins and in silico analysis of mutation impact on clinical phenotype and survival analysis of associations between mutations and hazard of clinical end-points. Fifteen SOD1 mutations were identified in 21.1% familial and 2.3% sporadic ALS cases. Their effects on SOD1 protein structure and functioning inferred from molecular modeling and in silico analyses correlate well with the clinical data. Molecular modeling results support the hypothesis that folding intermediates rather than mature SOD1 protein give rise to the source of cytotoxic conformations in ALS. Significant associations between type of mutation and clinical end-points were found.

The EcCLC antiporter embedded in lipidic liquid crystalline films - molecular dynamics simulations and electrochemical methods.

Lipidic-liquid crystalline nanostructures (lipidic cubic phases), which are biomimetic and stable in an excess of water, were used as a convenient environment to investigate the transport properties of the membrane antiporter E. coli CLC-1 (EcCLC). The chloride ion transfer by EcCLC was studied by all-atom molecular dynamics simulations combined with electrochemical methods at pH 7 and pH 5. The cubic phase film was used as the membrane between the chloride donor and receiving compartments and it was placed on the glassy carbon electrode and immersed in the chloride solution. Structural characterization of lipidic mesoscopic systems with and without the incorporation of EcCLC was performed using small-angle X-ray scattering. The EcCLC transported chloride ions more efficiently at more acidic pH, and the resistance of the film decreased at lower pH. 4,4-Diisothiocyanatostilbene-2,2'-disulfonic acid (DIDS) employed as an inhibitor of the protein was shown to decrease the transport efficiency upon hydrolysis to DADS at both pH 7 and pH 5. The molecular dynamics simulations, performed for the first time in lipidic cubic phases for EcCLC, allowed studying the collective movements of chloride ions which can help in elucidating the mechanism of transporting the ions by the EcCLC antiporter. The protein modified lipidic cubic phase film is a convenient and simple system for screening potential inhibitors of integral membrane proteins, as demonstrated by the example of the EcCLC antiporter. The use of lipidic cubic phases may also be important for the further development of new electrochemical sensors for membrane proteins and enzyme electrodes.

Discovery of thiazolidin-4-one analogue as selective GSK-3ß inhibitor through structure based virtual screening.

GSK-3ß directly phosphorylate tubulin binding site of tau protein, indicating its importance in tau aggregation and, therefore, in Alzheimer's disease pathology. New GSK-3ß inhibitors were identified using a structure-based screening, ADMET analysis. These studies revealed that ZINC09036109, ZINC72371723, ZINC72371725, and ZINC01373165 approached optimal ADMET properties along with good MM-GBSA dG binding. Protein kinase assays of these compounds against eight disease-relevant kinases were performed. During disease-relevant kinase profiling, ZINC09036109 ((E)-2-((3,4-dimethylphenyl)imino)-5-(3-methoxy-4-(naphthalen-2-ylmethoxy)benzyl)thiazolidin-4-one) emerged as a selective GSK-3ß inhibitor with more than 10-fold selectivity over other disease-relevant kinases. Molecular dynamics study of ZINC09036109 molecule revealed interactions with Ile62, Phe67, Val135, Leu188, Asp200 amino acid residues of the binding site of GSK-3ß, which were highly comparable to the co-crystallized molecule and hence validating comparative better activity of this compound compared to overall screened molecules.

Unexpected Reaction Products of Uracil and Its Methyl Derivatives with Acetic Anhydride and Methylene Chloride.

New acetyl derivatives of uracil, 6-methyluracil, and thymine were obtained in the course of an unconventional synthesis in methylene chloride. It was shown that products with the acetyloxymethyl fragment are formed according to a mechanism different from that for products with the acetyloxyethyl group. In particular, for uracil it was proven that the reaction with Ac2O, TEA, and CH2Cl2 leads to 1-acetyloxymethyluracil, where the N1 substituent is composed of the -CH2- fragment that originated from CH2Cl2 and the 1-acetyloxy moiety from Ac2O. The reaction of uracil with Ac2O, TEA, CH2Cl2, and DMAP leads to an acetyloxyethyl derivative in which the -CH2-CH2- fragment originates from TEA and the 1-acetyloxy moiety from Ac2O. A possible mechanism for the formation of new compounds was suggested and supported by the density functional theory/B3LYP quantum mechanical calculations. New compounds (39 in total, including seven deuterated) were fully characterized by nuclear magnetic resonance and high-resolution mass spectrometry techniques.

The Role of Cholesterol in Amyloidogenic Substrate Binding to the γ-Secretase Complex.

Alzheimer's disease is the most common progressive neurodegenerative disorder and is characterized by the presence of amyloid ß (Aß) plaques in the brain. The γ-secretase complex, which produces Aß, is an intramembrane-cleaving protease consisting of four membrane proteins. In this paper we investigated the amyloidogenic fragments of amyloid precursor protein (substrates Aß43 and Aß45, leading to less amyloidogenic Aß40 and more amyloidogenic Aß42, respectively) docked to the binding site of presenilin, the catalytic subunit of γ-secretase. In total, we performed 9 µs of all-atom molecular dynamics simulations of the whole γ-secretase complex with both substrates in low (10%) and high (50%) concentrations of cholesterol in the membrane. We found that, at the high cholesterol level, the Aß45 helix was statistically more flexible in the binding site of presenilin than Aß43. An increase in the cholesterol concentration was also correlated with a higher flexibility of the Aß45 helix, which suggests incompatibility between Aß45 and the binding site of presenilin potentiated by a high cholesterol level. However, at the C-terminal part of Aß45, the active site of presenilin was more compact in the case of a high cholesterol level, which could promote processing of this substrate. We also performed detailed mapping of the cholesterol binding sites at low and high cholesterol concentrations, which were independent of the typical cholesterol binding motifs.

Homology Modeling Using GPCRM Web Service.

Homology modeling methods are commonly used for quick and precise construction of a desired protein or its mutant using protein templates, which were determined by crystallography, cryo-EM, or NMR. Due to the increasing number of such structures, the obtained models are precise even in the case of small similarity between sequences of template and modeled proteins. The reason for that is a high evolutionary conservation in the structure regions responsible for keeping the function of proteins. This is also the case for G protein-coupled receptors (GPCRs), which constitute the largest family of membrane receptors with nearly 800 proteins. The GPCRM web service ( https://gpcrm.biomodellab.eu/ ) was set up for the nearly automatic generation of high-quality structures of modeled GPCRs. The three possible paths: "High similarity," "Quick path," and "Long path" allow the user to choose between a fast but less reliable path, up to more reliable but longer procedures. In the Advanced mode the service allows for user modifications including selection of template(s) and a manual adjustment of the sequence alignment.

GPCRsignal: webserver for analysis of the interface between G-protein-coupled receptors and their effector proteins by dynamics and mutations.

GPCRsignal (https://gpcrsignal.biomodellab.eu/) is a webserver devoted to signaling complexes of G-protein-coupled receptors (GPCRs). The recent improvement in cryo-electron microscopy resulted in the determination of a large number of high-resolution structures of GPCRs bound to their effector proteins: G proteins or arrestins. Analyzing the interfaces between receptor and an effector protein is of high importance since a selection of proper G protein or specific conformation of arrestin leads to changes of signaling that can significantly affect action of drugs. GPCRsignal provides a possibility of running molecular dynamics simulations of all currently available GPCR-effector protein complexes for curated structures: wild-type, with crystal/cryo-EM mutations, or with mutations introduced by the user. The simulations are performed in an implicit water-membrane environment, so they are rather fast. User can run several simulations to obtain statistically valid results. The simulations can be analyzed separately using dynamic FlarePlots for particular types of interactions. One can also compare groups of simulations in Interaction frequency analysis as HeatMaps and also in interaction frequency difference analysis as sticks, linking the interacting residues, of different color and size proportional to differences in contact frequencies.

Allosteric Modulation of the CB1 Cannabinoid Receptor by Cannabidiol-A Molecular Modeling Study of the N-Terminal Domain and the Allosteric-Orthosteric Coupling.

The CB1 cannabinoid receptor (CB1R) contains one of the longest N termini among class A G protein-coupled receptors. Mutagenesis studies suggest that the allosteric binding site of cannabidiol (CBD) involves residues from the N terminal domain. In order to study the allosteric binding of CBD to CB1R we modeled the whole N-terminus of this receptor using the replica exchange molecular dynamics with solute tempering (REST2) approach. Then, the obtained structures of CB1R with the N terminus were used for ligand docking. A natural cannabinoid receptor agonist, Δ9-THC, was docked to the orthosteric site and a negative allosteric modulator, CBD, to the allosteric site positioned between extracellular ends of helices TM1 and TM2. The molecular dynamics simulations were then performed for CB1R with ligands: (i) CBD together with THC, and (ii) THC-only. Analyses of the differences in the residue-residue interaction patterns between those two cases allowed us to elucidate the allosteric network responsible for the modulation of the CB1R by CBD. In addition, we identified the changes in the orthosteric binding mode of Δ9-THC, as well as the changes in its binding energy, caused by the CBD allosteric binding. We have also found that the presence of a complete N-terminal domain is essential for a stable binding of CBD in the allosteric site of CB1R as well as for the allosteric-orthosteric coupling mechanism.

Molecular Modeling of Histamine Receptors-Recent Advances in Drug Discovery.

The recent developments of fast reliable docking, virtual screening and other algorithms gave rise to discovery of many novel ligands of histamine receptors that could be used for treatment of allergic inflammatory disorders, central nervous system pathologies, pain, cancer and obesity. Furthermore, the pharmacological profiles of ligands clearly indicate that these receptors may be considered as targets not only for selective but also for multi-target drugs that could be used for treatment of complex disorders such as Alzheimer's disease. Therefore, analysis of protein-ligand recognition in the binding site of histamine receptors and also other molecular targets has become a valuable tool in drug design toolkit. This review covers the period 2014-2020 in the field of theoretical investigations of histamine receptors mostly based on molecular modeling as well as the experimental characterization of novel ligands of these receptors.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.